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We report a complete calculation of the quark and glue momenta and angular momenta in the pro- 
ton. These include the quark contributions from both the connected and disconnected insertions. 
The calculation is carried out on a 16 3 x 24 quenched lattice at /3 =6.0 and for Wilson fermions 
with K = 0.154,0.155, and 0.1555 which correspond to pion masses at 650, 538, and 478 MeV. 
The quark loops are calculated with Z\ noise and signal-to-noise is improved further with un- 
biased subtractions. The glue operator is comprised of gauge-field tensors constructed from the 
overlap operator. The u and d quark momentum/angular momentum fraction is 0.66(5)/0.72(5), 
the strange momentum/angular momentum fraction is 0.024(6)/0. 023(7), and that of the glue is 
0.3 1(6)/0. 25(8). The orbital angular momenta of the quarks are obtained from subtracting the 
angular momentum component from its corresponding spin. As a result, the quark orbital angular 
momentum constitutes 0.50(2) of the proton spin, with almost all it coming from the disconnected 
insertion. The quark spin carries a fraction 0.25(12) and glue carries a fraction 0.25(8) of the total 
proton spin. 
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1. Introduction 

Determining the contributions of quarks and gluons to the nucleon spin is one of the challeng- 
ing issues in QCD both experimentally and theoretically. Since the contribution from the quark 
spin is small (~ 25%) from deep inelastic scattering experiments, it is expected that the rest should 
come from glue spin and the orbital angular momenta of quarks and glue. 

Lattice calculations of the quark orbital angular momenta have been carried out for the con- 
nected insertions [jl[ ^, |3| |J and it was shown to be small [jl]| in the quenched calculation and near 
zero in dynamical fermion calculations [Q, ||, |[] due to the cancellation between those of the u and 
d quarks. Gluon helicity distribution AG(x)/G(x) from both COMPASS and STAR experiments is 
found to be close to zero [gj. Furthermore, it is argued based on analysis of single-spin asymmetry 
in unpolarized lepton scattering from a transversely polarized nucleon that the glue orbital angular 
momentum is absent [0], Thus, it appears that we have encountered a 'Dark Spin' scenario. 

In this work, we give a complete calculation of the quark and glue momenta and angular mo- 
menta. The quark contributions in both the connected and disconnected insertions are included. 
Combining with earlier work of the quark spin, we obtain the quark orbital angular momenta. We 
find that indeed the u and d quark orbital angular momenta largely cancel in the connected inser- 
tion. However, their contributions including the strange quark are large (50%) in the disconnected 
insertion due to the fact that the quark spin for each of the u, d, and s quarks in the disconnected 
insertion is large and negative. We have been able to obtain the glue momentum and angular mo- 
mentum for the first time, mainly because the overlap operator is used for the gauge-field tensor 
which is less noisy than that constructed from the gauge links. The lattice renormalization of the 
quark and glue energy-momentum operators are achieved through the momentum and angular mo- 
mentum sum rules. Due to these sum rules, the gravitomagnetic moment is zero as proven for 
composite systems from the light-cone Fock representation [^]. 

2. Formalism 

2.1 Angular Momenta and Momenta for Quarks and Gluons 

The angular momentum operator in QCD can be expressed as a gauge-invariant sum []|], 

jQCD=Jq+Jg=-Z^ + Lq + J g (2.1) 

where J q = ^L + L q and J g are the quark and gluon contributions respectively. £ is the quark spin 
operator, and L q is the quark orbital angular momentum operator. J q>g can be expressed in terms of 
the energy-momentum tensor operator through 

Jf' 8 = \ etjk J d 3 x(3™xi - 3Jf **) (2.2) 

Similarly, the quark and glue momentum operators are 

Pf'*= f d?x2T™, (2.3) 



Quark and glue momenta and angular momenta 



K.F. Liu 



where 



%v = X L V/[74 A +7* A -74 A- -Yi D 4 ]Vf 



(2.4) 



and 



«55f = i 2^ g m Gut 

k=\ 



(2.5) 



The matrix element of the energy-momentum tensor, SFm, between two nucleon states with mo- 
menta, p' and p, can be written as [Ist] (in Euclidean space), 



(p,s\&™\ p ', s 



u(p,s) 



Tl(q 2 )(Y4Pi + YiP4) ~ ^—T 2 (q 2 )(p40iaq a +Pi(y4aqa) 
Zm 



111 



-Ti(q 2 )q4qt 



u(p',s' 



(2.6) 



<i-x 



where, p = ^(p+p'), q^= p^—p'mh the mass of the nucleon, and u(p,s) is the nucleon spinor. 
By substituting Eq. (Elj) into Eqs. (2.2) and (2.3) at q 2 — > limit, one obtains 



h, g = ^[Tl(0) + T2(0)] q>g , 



(x) q , g = 7i(0) M . 



(2.7) 
(2.8) 



In this present study, there are two operators — energy momentum tensors for the quarks and 
glue. We can use the momentum and angular momentum sum rules to perform the renormalization 
on the lattice 



Z ? (a)ri(0) ? +Z g (a)7i(0) g = l, 

Z q {a)[T l (0) q + T 2 {0) q ]+Z g {a)[T l (0) g + T 2 {0) g ] = l. 



(2.9) 
(2.10) 



After the lattice renormalization, one can perform perturb ative mixing and matching to the ms 
scheme at /J. = 2 GeV. 

It is interesting to note that from these equations, i.e. Eqs. (2.5) and ( 2.1CJ ), one obtains the 
sum of the 72(0) 's for the quarks and glue to be zero, i.e. 



Z q (a)T 2 (0) q +Z g (a)T 2 (0) g =0. 



(2.11) 



This shows that the total anomalous gravitomagnetic moment of the nucleon vanishes. This has 
been proven by Brodsky et al. ^ for composite systems from the light-cone Fock representation 
and now it is shown as a consequence of momentum and angular momentum conservation. We will 
see how large the individual quark and glue contributions are. 

Since T\(q 2 ) and T 2 (q 2 ) can have different q 2 behaviors, we shall compute them separately at 
different q 2 values []2p, and then separately extrapolate them to q 2 — > for both the quark and glue 
contributions. 
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2.2 Glue Energy-Momentum Tensor Operator 

It is well-known that gauge operators from the link variables are quite noisy. We adopt the 
glue energy-momentum tensor operator where the field tensors are obtained from the overlap Dirac 
operator. The gauge-field tensor has been derived from the massless overlap operator D ov [P] 



Tr s [a MV D 0V (x,x)] = c T a 2 G^ v (x) + &(a i ) 



(2.12) 



where Tr s is the trace over spin, cj = 0.11157 for K = 0.19. With this construction, we expect 
the ultraviolet fluctuations to be suppressed, since D ov is exponentially local and D ov (x,x) involves 
gauge loops beginning and ending at x, which serves as smearing. 



3. Lattice Calculations and Numerical Parameters 

In order to obtain J qig = j [7\ (0) + ^(0)] and {x) q _ g = T\ (0) 9:g , we first calculate the three- 
point functions, GN^r 4i N(p,t 2 ',q,ti;p' ,to) and the two-point functions G^N{p,h)- Here p is the sink 
momentum, q is the momentum transfer, and p' is the source momentum, to, t\ and t 2 are the 
source, current insertion, and sink time respectively. We then take the following ratios between 
three-point and two-point functions, which involve the combinations of T\{q 2 ),T 2 {q 2 ), and Ti,{q 2 ) 



Tr 



r e , m G/v^ 4/ ^(^, ti\ —q, t\ ;//, to) 



Tr[r e G NN (p,t 2 )] 



Tr 



T e G NN (p',t 2 -ti+to) 



\ Tr[r e G NN (p,t 2 -ti+to)] 



Tr [T e G m (p, h )] Tr [r e G NN (p, t 2 )] 



\ Tr \T e G NN {p',h )J Tr \T e G NN {p' , t 2 ) 



fl>?0,f2>fl 



[aiTijq 2 ) + a 2 T 2 (q 2 ) + a^T^q 2 
AyjE p {E p + m)E p , (Ep* + m) 



-(3.1) 



where Y m = (— j)/2(l + ya)y m y^ is the spin polarized projection operator, and Y e = 1/2(1 + 74) is 
the unpolarized projection operator. For the special case with p = 0, one obtains the combination 

ofin + T^iq 2 ) 



Tr 



r m Gjv5yv(0,? 2 ;-<7,?i;/?' ,t ) 



Tr 



F e GNN{0,ti] 



Tr 



r e GAw(0,?2) 



Ti[T e G m (-q,h)] «»^2»« 



£ij m qj[Ti+T 2 }{q 2 ), (3.2) 



and the forward matrix element, which is obtained with the following ratio 



Tr 



r e G N ,3r 4jN (p,t 2 ;0,ti;p,to) 



\p\Tr[F e G NN (p,t 2 )} 



'1>'0.'2>'1 



(x). 



(3.3) 



Since the quark T\{q 2 ) and T 2 (q 2 ) for the connected insertion are shown to have quite different 
q 2 behavior ^ ||, Q], we will need to separately extrapolate T\{q 2 ) and T 2 {q 2 ) to q 2 — > 0. We 
combine several kinematics and with both the polarized and unpolarized three-point functions to 
extract T\(q 2 ) ,T 2 {q 2 ) , and Ti,{q 2 ) which appear as different combinations in the 3-point to 2-point 
ratios in Eqs. ( p.l[ ) and then extrapolate 7\ (q 2 ) and T 2 {q 2 ) in q 2 to obtain T\ (0) and T 2 {§). We can 
check this by comparing the extracted T\ (q 2 ) + T 2 (q 2 ) against the [7\ + T 2 ] (q 2 ) obtained directly 
from Eqs. (EO) and (3.3) at comparable q 2 . 
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The three-point functions for quarks have two topologically distinct contributions in the path- 
integral diagrams — one from connected (CI) and the other from disconnected insertions (DI) | JTo| , 
For DI, we sum over the current insertion time, t\, between the source and the sink 



11, 121 



time, i.e. from t\ = to + 1 to ?2 ~ 1 [jl3|» |14[ |15|, |16j] in order to increase the statistics. Similarly 



for the glue. Then the corresponding ratios at large time separation for Eqs. (3.1), (p.2[), and 



\> are 4 / Ep (E p +m)E pl (E pl+ m) Xh + const -> %««; [ T l + T 2\ {T)q,g Xt 2 + const., and (x) Xt 2 + const. 

respectively. 

We then extract the slopes and obtain T l (0),T 2 (0), \T\ + T 2 ] {q 2 ) q , g and {x) q>g in the DI the 
same way as was done for the CI. 

We use 500 gauge configurations on a 16 3 x 24 lattice generated with Wilson action at /3 =6.0 
in the quenched approximation. The values of the hopping parameter we use are k = 0.154,0.155, 
and 0.1555. The critical hopping parameter, K c = 0.1568 is obtained by a linear extrapolation 
to the zero pion mass Jl7|]. Using the nucleon mass to set the lattice spacing at a = 0.11 fm, 
the corresponding pion masses are 650(3), 538(4), and 478(4) MeV, and the nucleon masses are 
1291(9), 1159(11), and 1093(13) MeV, respectively. We use Dirichlet boundary condition in the 
present work. 

In the case of DI, the quark loop is evaluated separately. We compute it stochastically by 



using complex Z 2 noise vector Q18J. The number of noise vectors we use is 500 on each gauge 
configuration. Also for the case of quarks, we shall define two k:'s for the quark mass: K v for 
valence quarks, and fc sea for sea quarks. For the strange quark current, we have fixed K" sea = 0.154 
which is close to the strange quark mass as determined from the <p meson mass, and K v takes the 
values of 0.154,0.155, and 0.1555. For up and down quarks, we consider the cases with equal 
valence and sea quark masses, i.e. K" sea = K v = 0.154,0.155, and 0.1555. 

In the case of glue, we construct the energy-momentum tensor from the gauge field tensors 



which are calculated from the massless overlap operator, D ov H [l^] as discussed in Sec. 12. 



D ov (x,x) at all space-time points are estimated stochastically where we compute the color and spin 
indices exactly, but we perform space-time dilution by a separation of two sites on top of odd/even 
dilution. Therefore, the "taxi-driver distance" = 4 in our case. We use two Z\ noise sources, and 
take the average of them on each configuration. 

Due to the stochastic estimation, there is noise in addition to that from the gauge configura- 
tions. So, we adopt the following techniques to reduce the error: 

• We use discrete symmetries to discard the unwanted part of the current and the two-point 
functions when we correlate the current to the two-point functions to construct the three- 
point functions [20, ^TJ p. Combining 75-Hermiticity, parity and CH transformation, we 



determine whether the real or imaginary part of the current and the two-point functions will 
contribute to the three-point functions. 

We employ unbiased subtraction to the noise estimation of the quark loop to reduce the 
contributions from the off-diagonal matrix elements [^2|, [TL 15]- We use four subtraction 



terms (kD, k 2 D 2 , k 3 D 3 and K 4 D 4 ). No subtraction has been implemented for the estimation 
of the glue matrix elements. 
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We use multiple nucleon sources (in our case, 16) to increase the statistics. We correlate 
the nucleon propagators at different source locations with the already computed quark loops 



which results in significant reduction of errors [|5| |16| ]. 



The error analysis is performed by using the jackknife procedure. The correlation among 
different quantities are taken into account by constructing the corresponding covariance matrices. 
In order to extract various physical quantities, we use correlated least-^ 2 fits. 

4. Numerical Results 



We first present our results of the CI. In Fig. 1(a) we plot T\ u (q) + T 2iU (q ) and T\^{q ) + 



T2,d(q 2 ) as a function of — q 2 for the case of K = 0.1555, where T\(q 2 ) and T^ 2 ) are obtained 
from Eq. Q). We also plot [71 )B + T 2)U ] (q 2 ) and [T u + T 2>d ] (q 2 ) directly from Eq. (Q at slightly 
different — q 2 . We see that the latter basically agrees with the error band of the former with a dipole 
fit in — q 2 within 2 sigmas. This is a cross check of our procedure of extracting T\{q 2 ) and T 2 {q 2 ) 
from Eq. (3.1) which involves 3 to 4 equations of the combinations of T\{q ),T2(q 2 ), and T?,{q 2 ) 
at different q 2 . We also show, in Fig. |l(b)| , T\ M {q 2 ) + T\j(q 2 ) and T 2>u (q^) + 72,c/(<? 2 ) and their 



error bands. Also plotted is 71«(0) + 71^(0) from Eq. (3.3). We see that its error is smaller than 
that from the separately extrapolated 71„(0) and 71^(0). Thus we shall use 71 )M (0) + 71,^(0) from 
Eq. (3.3) and combine with T 2>u (q ) + T 2i d(q 2 ) to get the angular momentum for the CI. 



T,+T 2 , kappa=0.1555 (CI) 



T u (q 2 )+T d (q 2 ),kappa=0.1555(CI) 



0.8 
0.6 
0.4 
0.2 

- 
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Ti 4 (q 2 )tT 2 ;d(q -£ 
lJi,u +T 2,u](qp) ' — 
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-q 2 
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Figure 1: Results from K = 0.1555: (a) The sum of T\{q ) and T%(q ) extracted from Eq. ( J3 . 1[ > with error 



bands from the dipole fit is compared to [71 + T 2 ] (q 2 ) from Eq. (3.2) with slightly different — q 2 for u and 
d quarks in the CI. (b) The sum of u and d quark contributions for T\ (q 2 ) and T2(q 2 ). The blue asterisk at 
-q 2 = is 71,„(0) + r u/ (0) from Eq. (O). 



For the DI with K v = KT sea = 0.1555 at^ 2 a 2 = —0.144, we show in Fig. |2(a)| the ratio of Eq. (| 
with t\ summed between to + 1 and t 2 — l and plotted against the sink time t 2 so that the slope is 
£ijmqj [71 + T2] (q 2 ). We fit the slope from t 2 = 8 where the ratio is dominated by the nucleon to t 2 = 
12. We plot [71 + T 2 ] (q 2 ) so obtained in Fig. [2(b) and compare them to T\(q 2 ) + T2{q 2 ) extracted 
from 5 to 6 combinations of ci\T\(q 2 ) +fl , 27 , 2(g 2 ) +a^,Tj,{q 2 ). We see that they are consistent with 
each other within errors. The error bands are from the dipole fits of T\{q 2 ) and T2(q 2 ). 71(0) 
(red square) is from the forward matrix element which has smaller error than the — q 2 extrapolated 
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7i(0). We shall combine it with the extrapolated ^(0) for the angular momentum 7 in the DI. 
Finally, we perform a linear chiral extrapolation to obtain 7\(0) + ^(0) for the u/d quark at the 
chiral limit. This is shown in Fig. 2(c) . For the strange, we fix the quark loop at K" sea = 0.154 and 
extrapolate K v to the chiral limit. 

We perform similar analysis for the glue momentum and angular momentum. They are plotted 
in Figs. |(a| f(bj| and gcj, 



With all the quark and glue momenta and angular momenta, we carry out lattice renormal- 
ization through the sum rules in Eqs. (2.9) and (2.10). We obtain Z q (a) = 1.05 and Z g (a) = 1.05. 



This shows that both the lattice operators, particularly the glue energy-momentum tensor from the 
overlap operator, are 'natural' and close to the continuum. In Table [j], we list the quark momentum 
fractions (x) = T\ (0) in the CI (u and d) and the DI (u/d and s) and that of the glue. We also list 
the corresponding ^(0) and angular momentum 27 = T\ (0) + ^(0). 

Table 1: Lattice renormalized values with renormalization constants Z„=Z„=1.05. 





CI(m) 


Cl(d) 


Cl(u + d) 


Dl(u/d) 


Dl(s) 


Glue 


(x) 


0.428(40) 


0.156(20) 


0.586(45) 


0.038(7) 


0.024(6) 


0.313(56) 


r 2 (o) 


0.297(112) 


-0.228(80) 


0.064(22) 


-0.002(2) 


-0.001(3) 


-0.059(52) 


u 


0.726(118) 


-0.072(82) 


0.651(50) 


0.036(7) 


0.023(7) 


0.254(76) 


gA 


0.91(11) 


-0.30(12) 


0.61(8) 


-0.12(1) 


-0.12(1) 


— 


2L 


-0.18(16) 


0.23(15) 


0.04(9) 


0.16(1) 


0.14(1) 


— 



We see from Table [j] that the strange momentum fraction (x) s = 0.024(6) is in the range of 
uncertainty of (x) s from the CTEQ fitting of the parton distribution function from experiments 
which is 0.018 < (x) s < 0.040 Q23Q. The glue momentum fraction of 0.313(56) is smaller than, 
say, the CTEQ4M fit of 0.42 at Q = 1.6 GeV [||]. We expect the glue momentum fraction to be 
larger than the present result when dynamical configurations with light fermions are used in the 
calculation. From Figs. 1(b) and 3(b) and Table |], we find that the central values of ^(0) for 
u/d and s in the DI are small and consistent with zero. T2 M (0) + T2 j (0) in CI is positive, while 
?2,g(0) is negative. With the renormalization constants Z q (q) and Z g (a) fitted to be both positive 
(and near unity), they cancel giving a null anomalous gravitomagnetic moment. The fact that their 
unrenormalized magnitudes are almost identical is consistent with the finding that Z q (q) and Z g (a) 
are the same in size within errors. 

In analogy to 7^(0) which is known as the anomalous magnetic moment of the nucleon, ^(0) 
is termed the anomalous gravitomagnetic moment and has been shown by Brodsky et al. to vanish 
for composite systems |0j. We have now verified that this is the consequence of the momentum 
and angular momentum sum rules. 

The flavor-singlet g A , which is the quark spin contribution to the nucleon, has been calculated 



before on the same lattice [14]. We can subtract it from the total angular momentum fraction 27 to 
obtain the orbital angular momentum fraction 2L for the quarks. As we see in Table [j], the orbital 
angular momentum fractions 2L for the u and d quarks in the CI have different signs and they add 
up close to zero, 0.04(9). This is the same pattern seen with dynamical fermions configurations 
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light quark, [T 1+ T 2 ](q^) (Dl) 
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light quark(DI), chiral extrapolation 




-0.01 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 

m a 

(c) 



Figure 2: (a) The ratio of Eq. (3.2) withfi summed between to + 1 andf2~ 1 as a function of the sink time ti. 
The slope is fitted to obtain Sijmlj [Ti + T2] {q 2 ) q - (b) Separately extracted T\ (q 2 ) and Tiiq 2 ) are compared 
with [T\ + T2] (q 2 )- T\ (0) (red square) is from the forward matrix element, (c) Chiral extrapolation of T\ (0) 
and 7^(0) f° r tne u /d quark. They are not renormalized. 
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Figure 3: The same as in Fig. ||for the glue. 



Quark and glue momenta and angular momenta K.F. Liu 



with light quarks [Q, ^ Qj. The large 2L for the u/d and s quarks in the DI is due to the fact that 
g° A in the DI is large and negative, i.e. —0.12(1) for each of the three flavors. All together, the 
quark orbital angular momentum constitutes a fraction 0.50(2) of the proton spin. The majority 
of it comes from the DI. The quark spin fraction of the nucleon spin is 0.25(12) and glue angular 
momentum contributes a fraction 0.25(8). 

5. Summary 

In summary, we have carried out a complete calculation of the quark and glue momentum 
and angular momentum in the nucleon on a quenched 16 3 x 24 lattice with three quark masses. 
The calculation includes both the connected insertion (CI) and disconnected insertion (DI) of the 
three-point functions for the quark energy-momentum tensor. We used complex Z2 (or Z4) noise 
to estimate the quark loops in the DI and the gauge-field tensor from the overlap operator in the 
glue energy-momentum tensor. We find that we can obtain reasonable signal for the glue operator 
constructed from the overlap Dirac operator. After chiral extrapolation, we use the momentum and 
angular momentum sum rules to determine the lattice renormalization constants which turn out 
to have a 'natural' size close to unity, Z q {a) = 1.05 and Z g (a) = 1.05. The lattice renormalized 
momentum fractions for the quarks are 0.586(45) for the CI and 0.100(12) for the DI. The glue 
momentum fraction is 0.313(56). We have shown that the anomalous gravitomagnetic moment 



Eq. ( [2.1 1| ) vanishes due to the momentum and angular momentum conservation. 

After subtracting from the angular momentum 2/ the quark spin (g^) from a previous calcula- 



tion on the same lattice [|14j], we obtain the orbital angular fraction 2L. In the CI, we find that the u 
quark contribution is negative, while the d quark contribution is positive. The sum is small, 0.04(9). 
This behavior is about the same as observed in dynamical calculation with light quarks [Q, [| |J. 
The majority of the quark orbital angular momentum comes from the DI, because the quark spin 
from the DI is large and negative for each of the three flavors. In the end, we find the quark orbital 
angular momentum, the quark spin, and glue angular momentum fractions of the nucleon spin are 
0.50(2), 0.25(12), and 0.25(8) respectively. 

We are in the process of calculating the perturbative matching to the ms scheme with mixing 
so that we can quote our results in ms scheme at 2 GeV. 
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